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We have developed an all particle Fock-space relativistic coupled-cluster method for two-valence 
atomic systems. We then describe a scheme to employ the coupled-cluster wave function to calcu- 
late atomic properties. Based on these developments we calculate the excitation energies, magnetic 
hyperfine structure constants and electric dipole matrix elements of Sr, Ba and Yb. Further more, 
we calculate the electric quadrupole HFS constants and the electric dipole matrix elements of Sr + , 
Ba + and Yb + . For these we use the one- valence coupled-cluster wave functions obtained as an inter- 
mediate in the two-valence calculations. We also calculate the magnetic dipole hyperfine structure 
constants of Yb + . 
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I. INTRODUCTION 

Coupled-cluster theory, first developed in nuclear 
many body physics [TJ [2] , is considered one of the best 
many body theory. In recent times, it has been used with 
great success in nuclear [3] , atomic [H [5] , molecular [5] 
and condensed matter [7] calculations. A recent review 
[5] provides a detailed overview of the theory and varia- 
tions suitable to different classes of many-body systems. 
An earlier review provides an overview on the applica- 
tion of coupled-cluster theory to various areas of physics 
[9]. In atoms it is equivalent to incorporating electron 
correlation effects to all order. It has been used exten- 
sively in precision atomic structure and properties cal- 
culations. These include atomic electric dipole moments 
[H EPJ ; parity nonconservation , hyperfine structure 
constants [S] I12j and electromagnetic transition proper- 
ties (nun]. 

In this paper we report the development and results of 
relativistic coupled-cluster atomic calculations for two- 
valence atoms. For this we employ the Fock-space open- 
shell CCT [TSJ [TB] , which is also referred as valence uni- 
versal. Based on which the two-valence CC wave oper- 
ators are calculated via the closed and one- valence wave 
operators. The necessary developments of these inter- 
mediate stages of calculations were reported in our pre- 
vious works [T71 [IS]- We emphasize that in ref. [TSJ 
we proposed a new scheme to calculate properties with 
CC wave functions to all order. In the present work 
we implement two-valence Fock-space relativistic CCT 
with an all particle valence space. A similar approach 
was adopted in a previous work on relativistic coupled- 
cluster calculations of two- valence systems [19]. This is 
general enough for the precise wave function and prop- 
erties calculations of the low-lying levels of two-valence 
systems like the alkaline-earth metal atoms, Yb and Hg. 
In these systems, the low- lying levels arise from the ns 2 , 
ns(n — l)d and nsnp configurations. We show selecting a 
model space consisting of these configurations is incom- 
plete but quasi-complete. Advantage of quasi-complete 
model space is, it has all the virtues of a complete model 
space but one can circumvent the divergence associated 



with intruder states J20J in open-shell CCT. The calcula- 
tions presented in this work are based on coupled-cluster 
singles and doubles (CCSD) approximation. It was ini- 
tially formulated for molecular calculations [H] and used 
in atomic structure calculations to study the excitation 
energies of Li |22] . Later, the relativistic version was im- 
plemented to calculate structure and properties of high 
Z atoms and ions [5] 133] 121]. 

In the present work we apply the method we have de- 
veloped to calculate the wave functions of alkaline-Earth 
atoms Sr and Ba, and lanthanide atom Yb. All of these 
atoms are candidates of extremely precise experiments 
either for application oriented investigations or to probe 
fundamental laws of nature. Atomic Sr, which was re- 
cently cooled to quantum degeneracy [25j . is a strong 
contender of future optical clocks [26] [27] . Experiments 
on Bosc-Einstein statistics violations have used Ba as the 
target atom |28| and it is an ideal proxy, both for exper- 
imental [21] and theoretical [30J studies, of atomic Ra. 
An atom with large parity and time reversal violation 
effects [31) . and promising candidate for future exper- 
iments. Recently, parity violation was detected in Yb 
[32) and ongoing experiments could lead to unambiguous 
detection of nuclear anapole [33 . An exotic parity violat- 
ing nuclear moment, which can be detected only through 
atomic experiments. There are also proposals to mea- 
sure atomic electric dipole moment, a signature of parity 
and time violations, with novel techniques [34] [35] . It 
must be mentioned that several isotopes of Yb has been 
cooled to degeneracy |36j and could be employed in fu- 
ture precision measurements. Further more, atomic Yb 
in an optical lattice is a candidate of frequency standard 

[33. 

As mentioned before, to obtain the two-valence wave 
operator, we compute the closed-shell and one-valence 
wave operators in the intermediate steps. We take ad- 
vantage of this and use the one-valence wave operator 
to compute the hyper fine constants of Sr + , Ba + and 
Yb + . For the first two ions, we reported the magnetic 
hyperfine structure constants (HFS) in our previous pa- 
per [18] , so we compute only the electric quadrupole HFS 
constants. Whereas for Yb + we compute both the mag- 
netic dipole and electric quadrupole HFS constants. In 
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addition, we also compute the electric dipole transitions 
matrix elements. Like the neutral atoms, all the ions 
are under experimental investigations for various preci- 
sion measurements. For example, a single trapped 87 Sr + 
is a suitable frequency standard [38] . There are similar 
experiments with Yb + |39j as an alternative frequency 
standard. And it is one of the frequency standards in 
the laboratory measurement of temporal variation of fine 
structure constant [ID]. These are application oriented 
precision experiments. The other fascinating prospect 
is the observation of parity nonconservation in a single 
137 Ba+ |H]. 

The paper is divided into seven sections. In the next 
section, that is Section. II, we give a brief description of 
many-body perturbation theory (MBPT) for two- valence 
systems. It provides the minimal description of concepts 
pertinent to development of two-valence coupled-cluster 
theory. Section. Ill is a short writeup on closed-shell 
and one-valence CCT followed by derivation of the two- 
valence CCT in some detail. These are in the context of 
complete model space. Incomplete model space CCT is 
explained in Section. IV and atomic Yb is discussed as an 
example. Calculation of properties, HFS constants and 
electric dipole transition, with CC wave functions is the 
topic of Section. V. In Section. VI, the important details 
of implementing two- valence CCT is explained, however, 
with emphasis on physics. Finally, results and discus- 
sions are reported in Section. VII. In the paper, all the 
calculations and mathematical expressions are in atomic 
units (e = % = m e = 47reo = 1)- 

II. MBPT FOR TWO- VALENCE ATOMS 

Relativistic effects are the key to obtain accurate re- 
sults in the structure and properties calculations of high 
Z atoms with za ~ 1. The Dirac-Coulomb Hamiltonian 
is an approximate but an appropriate Hamiltonian 
to describe the properties of such atoms. For an atom 
with N electrons 

- 1 

h dg = [«* Pi + (P - i)c 2 - v N (n)] + E — ' 

i—1 i<j % ^ 

where a.^ and (3 are the Dirac matrices, p is the linear 
momentum, Vjv(r) is the nuclear Coulomb potential and 
last term is the electron-electron Coulomb interactions. 
It satisfies, in the case of two-valence atoms, the eigen 
value equation 

H BC \V VW ) =E va \¥ w ), (2) 

where indexes v and w represent the valence orbitals, 
l^mo) is the exact wave function and E vw is the exact 
energy of the two-valence atomic system. In MBPT, 
the total Hamiltonian, in Eq.Q, is separated into two 

parts: H = Y,ii ca i ' Pi + (« ~ l ) c2 ~ v NTi + w(n)], 
the unperturbed or exactly solvable part, and V — 



Si<j T~ ~ Ei u ( r i)j tne perturbation, referred as the 
residual Coulomb interaction. The unperturbed eigen 
functions \Q VW ) are the solution of the Dirac-Fock equa- 
tion 

H \$ VW )=EW\$ VW ), (3) 

here, | *f> vw ) are the antysymmetrised many-electron wave 
functions. Formally, in operator notations, these are 
generated from the closed-shell reference state |$o) as 

\$ vw ) = <4at|$o)- And, the eigen value, Evw, is the sum 
of the single electron energies. These are the basic start- 
ing points common to atomic MBPT and coupled-cluster 
theory (CCT). The two theories share a common thread 
till the generalized Bloch equation discussed in the 
next section, but is significantly different from there on. 

A. Generalized Bloch equation 

In this section, and others as well, we provide the ba- 
sic equations and necessary definitions essential to a lucid 
description of the method we have developed and used. 
For detailed descriptions appropriate references are pro- 
vided. In MBPT, the Hilbcrt space of the eigen functions 
is separated into two sub-manifolds: model space 
(P), comprise of the eigen functions which are the best 
approximation to the exact eigen functions of interest 
and remaining spans the orthogonal space (Q). In the 
single reference theory, the exact state \^ vw ) and model 
state \$ vw ) are related as 

\y vw } = n\$ vw ), (4) 

where, is the wave operator and is the solution of the 
generalized Bloch equation 

[n, h ]p = (vn - npvn)p. (5) 

Detailed exposition of the equation and relevant deriva- 
tions are given ref. j42j [43] • Here, intermediate normal- 
ization 

\$ vw ) = PSl\$ vw ), (6) 

is a necessary condition to obtain the generalized Bloch 
equation. In singles and doubles approximation, often 
used and well tested method in atomic calculations, the 
wave operator is 

il = 1 + x p a a) p a a + x p v ala v + \^x v * b a\a\a h a a + x p a q v ala\a v a a 

2%vw a p a q a w a v 7 (7) 

where ab ■ • • (pq ■ ■ •) denote core (virtual) orbitals and x.'.\ 
are the excitation amplitudes. This definition is crucial 
to our later discussions on the Fock space coupled-cluster 
in complete model space (CMS). Unlike the close-shell 
atoms, the model wave functions are not known in the 
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case of open-shell systems. These are obtained by diago- 
nalizing the effective Hamiltonian (iJ e ff ) matrix, calcu- 
lated within the P sub-manifold. Once the model wave 
function is obtained, the exact energy is the expectation 
value of -f/eff , as it satisfies the cigen value equation 



eff | ^VW / 



E |(E> ) 



(8) 



The effective Hamiltonian, in Eq.([8]), is expressed as 
H cS = PH P + PVSIP, 



(9) 



where Hq and V, as defined earlier, are the zeroth-order 
Hamiltonian and the residual Coulomb interaction re- 
spectively. The first term in Eq.Q is the leading contri- 
bution, E§w, to the total energy E vw . And the second 
term, with the wave operator f2, is the correction to 
referred as correlation energy. 



B. First- and second-order effective Hamiltonians 

From Eq.([8]), the first-order correction to energy is the 
expectation value of the first-order effective Hamiltonian 



o-U) _ 

— 



< eff = PVP = P(V 



Vi + V 2 )P. 



where, Vq is the contribution from the close-shell part and 
represented by the closed diagrams with no free lines at 
the vertexes. We exclude this term while calculating the 
excitation energies as it is common to all the diagonal 
elements of the -ff c ff matrix. It effectively shifts all the 
energy levels equally and does not account for the energy 
level splitting. The one- and two-body terms, V\ and 
Va, have contributions from open-shell part only. The 
contributing diagrams are the closed diagrams with free 
valence lines at the vertixes. The one-body term, Vi, 
also contributes to the diagonal elements only and hence 
doe s no t contribute to the energy level splitting. From 



Eq.(10l, H^g is reduced to the form 



H$ = PV 2 P. 



(11) 



This term contributes through a closed diagram with one 
pair of valence lines at each vertex shown in Fig. [2ja). 
From Eq.Q, the second-order effective Hamiltonian 



(a) 





(c) 



(d) 



FIG. 1: One-body diagrams arising from the second-order 
effective Hamiltonian H\* . 



h£> = pvnWp = p(Vi + y 3 )(fli + ^2 ) p , (i 2 ) 

where the superscript in the wave operator represents 
the order of the perturbation. Contributing diagrams 
are closed diagrams with valence orbitals as free lines 
at the vertexes. Detailed description of the relativistic 
second MBPT of two- valence systems is given ref. |44) . 
The terms involving V\ are zero if Dirac-Fock orbitals are 
used in the calculations. The expression of H cS is 



H. 



(2) 
eff 



pv 2 n i 2 1) p. 



(13) 



where {A ■ ■ ■ B} represents contraction between two op 



(2) J 

erators A and B. The diagrams of Hig , in Eq.(13), are 



(b) 



.0 



(e) 



(f) 




(10) FIG. 2: The two-body diagram (a), arises from the first- 



order effective Hamiltonian H^g . The remaining two-body 
diagrams, from (b) — (g), contribute to the second-order effec- 



(2) 



tive Hamiltonian H, 



separated into two categories. The diagrams with one 
pair of free lines as the valence orbitals, shown in Fig. [T] 
constitute the one-body effective operator. And the dia- 
grams with two pair of free lines as the valence orbitals, 
shown in Fig. [2j forms the two-body effective opera- 
tor. It must be mentioned that, a previous work reported 
the third ordered relativistic MBPT calculations of two- 
valence systems beryllium and magnesium iso-electronic 
sequences [4"5] . 



C. i^ ff matrix elements with jj coupled states 

In our scheme of calculations, we first evaluate the di- 
agrams arising from Eqs.(ll) and (13), Figs. [I] and (J2J, 
using uncoupled states. And we then store these as the 
effective one- and two-body operators. Later we use these 
effective operators to generate the matrix elements with 
respect to the jj coupled states. For two non- equivalent 
electrons the jj coupled antysymmetrised state may be 
expressed, in terms of the total angular momentum J 
state, as 



\{lvjvm v j w j w m w }JM) 



)JM) 



+ (-1) 



\(jwj w m w j v j v m v )JM} 



(14) 
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where, j v and j w are the total angular momenta of the 
single electron states \<fi v ) and \<p w ) respectively, 7^ and 
7u, are additional quantum numbers to specify the states 
uniquely. And m v and m w are the corresponding mag- 
netic quantum numbers. Similarly, J and M are the 
total angular momentum of the coupled state and mag- 
netic quantum number respectively. The matrix element 
of a two-body operator then consists of four terms, two 
direct and two exchange, with the normalization factor 
1/2. 

To evaluate the two-body matrix element, for exam- 
ple, the coulomb interaction shown Fig. [2^i. The direct 
matrix element is of the form 

((jvjvm v j w j w m w )JM\ — \(j x j x m x jyjyrriy)J'M') = 

T\1 



^2(-iy* +j ™+ J+k S(J, J')S(M, M') 



jx jx & 
jy jw J 



x(7^l|C' £ ||7 :c J a: )(7^l|C' c ||7 y J a > x R k 



(15) 



Where x and y represent valence orbitals, R k is the ra- 
dial integral and C k is the spherical tensor operator. 
For matrix elements, in Eq. (15), to be non-zero the 



states should have the same parity and J. The relation 
in Eq.( 15 ) holds true for the matrix elements of the other 
two-body diagrams Fig. [2^b-g) . In this case the multi- 
pole k and the radial integral arise from the combination 
of two orders of residual Coulomb interactions. 

Similarly, the matrix element of the one-body operator 
of rank k 7 with respect to the jj coupled state is 



{{lv3vl W J W )JM\F k {l)\( lxJxl yjy)j'M') = 

J k J' \ J j v j x k 
■M M )\ J' J j w 



(7^l|f fc H7xJx). (16) 



This is a very general expression and applicable to one- 
body operator of any rank k. In our calculations, how- 
ever, we use k = as the one-body effective operator is 
scalar. 



III. FOCK SPACE CCT: COMPLETE MODEL 
SPACE 

A model space is complete, if it consists of all the con- 
figurations formed by accommodating the valence elec- 
trons among the valence shells in all possible combina- 
tions. A remarkable consequence of choosing CMS in 
Fock-space coupled-cluster is that, the excitation opera- 
tors x'.W are common to all the determinants in the model 
space. Further more, x"\ uniquely separates into internal 
and external sectors. The external excitations contribute 
to O and projects a model function to the complemen- 
tary space. Whereas, internal excitations connect one 
model function to another model function and occurs in 
the definition of H e g. As we shall explore later, in the 



context of incomplete model space (IMS), such a neat 
separation is specific to CMS and another class of model 
space referred as quasi complete 03] • Validity of linked 
cluster theorem, one basic condition for any legitimate 
many-body theory, is assured in CMS. 
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FIG. 3: Low-lying energy levels of atomic Yb. 



The CMS, though endowed with several desirable 
properties, has one serious short coming for systems with 
two or more valence electrons. It inevitably encounters 
intruder states |20] and the outcome is severe convergence 
problems. This is the manifestation of model states with 
high energies that lies within the energy domain of the 
orthogonal space. In other words, in CMS when all possi- 
ble configurations are considered, the model and comple- 
mentary space are no longer energetically well separated. 
The occurrence of vanishing energy denominators is then 
a distinct possibility. Indeed, we invariably encounter 
intruder states in all our two-valence calculations with 
CMS. Its presence is the rule rather than the exception. 

For a better description of the CMS and intruder states 
let us consider a specific example, the low-lying levels 
of Yb atom. Configurations and terms of the ground 
and first few excited states important in precision spec- 
troscopy are 6s 2 ( 1 Sq), 6s6p ( 3 Pj), 5d6s ( 3 -Dj) and 
6s6p ( 1 -Pi). The 6s, 6p and 5d are then the obvious 
choice of valence shells. CMS of the system then consists 
of the configurations: 6s 2 , 6s6p, 5g?6s, 5d6p, 6p 2 and 5d 2 
and all the other configurations are in the complementary 
space. As shown in Fig. [3j the levels from the orthogonal 
space 6p7p { 3 Pj), 6s7s ( and 6s7s ( 1 S'o) lie within the 
model space. With several orthogonal functions within 
the energy domain of model functions, CMS based CCT 
calculations of Yb are likely to face with intruder state 
related divergences. Indeed, we do encounter divergences 
which, on careful analysis, can be attributed to the in- 
truder states. 
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In this work, we proceed to the relativistic two-valence 
coupled-cluster theory via the closed-shell [17] and one- 
valence coupled-cluster [TS] theories reported in our pre- 
vious works. We implement this within the framework 
of Fock-space or valence universal CCT 07]. The 
theory can be extended to systems with both particles 
and holes, however, for our present study an all parti- 
cle implementation is sufficient. Accordingly, the valence 
electrons are treated as particles [H] and each sector- 
closed-shell, one- and two-valence-are separate Hilbert 
spaces. Technical advantage of Fock-space CCT with 
CMS is the sector wise clean separation of cluster op- 
erators [IH]- However, the Hilbert spaces of two- valence 
subsumes the one-valence after a direct product with a 
spectator valence state and similarly, the closed-shell af- 
ter direct product with two-valence states. The state 
universal |50| is another flavour of open-shell CCT, where 
the wave operator is calculates in a single Hilbert space 
consisting of all the valence electrons. The wave operator 
is then state dependent and there is a lack of generality 
in the cluster equations. Further more, it requires com- 
plicated book keeping. These reasons have motivated us 
to choose the Fock-space CCT. Detailed discussions on 
Fock-space CCT and subtle issues related to the choice 
of model spaces are given in the review of Lindgren and 
Mukherjee [51]. 



Si 1 ' 



FIG. 4: Representation of the closed-shell and open-shell clus- 
ter operators 



A. Closed-shell and one-valence CCT 

Coupled-cluster theory is a non-perturbative many- 
body theory. It is equivalent to selecting linked terms 
in the Bloch equation, Eq.([5]), to all orders and combin- 
ing terms of same level of excitation (LOE). The elegance 
and perhaps, all the attendant difficulties of CCT is the 
exponential nature of the wave operator. We provide 
a brief reprise of the closed-shell [17] and one-valence 
CCT [TS], these form the initial steps of the two- valence 
coupled-cluster theory. Indeed, a large fraction of the 
cluster amplitudes in the two-valence theory arise from 
the close-shell. The exact atomic wave function of a one- 
valence system in the coupled-cluster theory is 



I*,: 



e T (l + S)\<S> v ) 



(17) 



where |$„) is the one- valence Dirac-Fock reference state, 
and T and S are the closed- and open-shell cluster op- 
erators respectively. As evident from the equation, all 
the higher order terms (non-linear) of T are calculated 



but S restricted to linear terms only. The later is on ac- 
count of the single valence electron. The diagrammatic 
representations of these operators are shown in Fig. [4j 

The closed-shell operator T in the coupled-cluster sin- 
gles doubles (CCSD) [5T] approximation is 



T = T 1 +T 2 , 



(18) 



where T\ and T 2 are the single and double excitation op- 
erator respectively. The closed-shell exact state in CCT 
is 



|*o) =e i |$ ) 



(19) 



and the cluster amplitudes are solutions of the nonlinear 
coupled equations 



ff N |$o) = 0, 



(*S|Fn^o>=0, 



(20) 
(21) 



where = e Hn e is the similarity transformed 
or dressed Hamiltonian. |$o) ; the Dirac-Fock reference 
state for the closed shell part, is the eigen value of the 
central potential Hamiltonian Hq. And and \^t) 
are respectively, the singly and doubly excited determi- 
nants. For details of the derivation, readers are referred 
to ref. [T7J. The open-shell cluster operator S is 



(22) 



■ And these are 



where and S< 2 > are the one-valence and two- valence 
cluster operators respectively. Similar to T, the open- 
shell one- valence cluster operator S^ 1 ', in CCSD approx- 
imation, is of the form — s[ 
solutions of the coupled linear equations 

($l\H N +{H^SW}\$ v ) = E* tt (<S>v\s{ 1} \$ v ), (23) 



In these equations E% tt is the attachment energy of an 
electron to the v shell. It is defined as 



ET = E„ 



(25) 



where E v = ($ V \H N + {H N S^}\^ V ) and E = 
($o\H\$o) these are the exact energy of \^ v ) and \^f ) 
respectively. The excited determinants, and |$^), 
are obtained by exciting an electron from valence orbitals 
to the virtuals. For detail description, of the derivation 
and interpretations, one may see r ef |18j. Nonzero renor- 
malization, right hand side in Eq. ( 23p4 ), is the predom- 
inant departure of open-shell CC from closed-shell CC. 
In the language of many-body diagrams, folded diagrams 
embody the renormalization terms. These, the folded 
diagrams, are topologically very different from the dia- 
grams of H or HnS^\ To illustrate the difference folded 
diagrams from the two-valence CC are shown in Fig. [5] 
Strictly speaking, the distortion in these diagrams are 
introduced to obtain correct energy denominators with 
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FIG. 5: Folded diagrams from the renormalization term 
in the generalized Bloch equation of two- valence systems. 
In two-valence coupled-cluster theory these diagrams arise 
from (a) E^(^ W \S^\^ VW ), (b)^2(*S8,| 1 S^ 1) |*„„> and 
(c)E^(& v l\Si 2) \* vw ). 



On the right hand side we have used the relations 
($ VW \S\$ VW ) = and ($> V w\S 2 \$> vw ) = 0, as the oper- 
ation of 5 on the state l^,™) transforms it to an excited 
determinant orthogonal to ($«u)|. 

To obtain the two- valence cluster equations, we project 
Eq. ( 28 ) on the doubly excited determinants 



AEZ 1 



1 + 5+ ^5 



1 4+ 



(30) 



the diagrammatic evaluation in MBPT. CC being non- 
perturbative there is no reason to be concerned about 
correct denominators. However, we retain the nomencla- 
ture and structure, in the diagrammatic analysis of CC 
equations, to identify the diagrams uniquely. 



Two-valence CCT 



In the two-valence sector, the cluster operator — 
S2 , a natural outcome of treating the single valence 
excitations as one-valence problem in Fock-space CCT. 
The exact two-valence state in CCT is 



+ s[ 1) + 1 -s{ 1) \si 1) + s : 



(2) 



(26) 



Here, for the valence part we have used exp(5) = 1 



S[ 1] + (1/2)S[ 1V 



S. 



(i) 



s. 



(2) 



Notice that, though 



(1/2)5} does not contribute to the one- valence CC 
equations, it does contribute to the two- valence CC equa- 
tions. Using this in Eq.(|2j) and projecting on e~ T , we get 



H 



1 + S[ r) + ^5} 



1 



5 



(i) 



(2) 



+5} 



(i) 



1 Mir 
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(i) 



5 



(2) 



5 



(27) 



Here after for simplicity of representation we use 
exp(5) = l + 5+(l/2)5 2 with the definition, restricted to 

two-valence sector only, 5 2 = 5^ . Using the normal- 
ordered form the Hamiltonian, H = 
write 

1 



E®£, we can 



1 



1 + 5+ -S' 



= AE° 



1 + 5 



S 2 



1$,- 



(28) 



where A£™ rr 



En, 



E^, is the correlation energy 
of the two- valence atoms and as defined earlier, Eq.(22|, 
5 = S^ 1 ' + S^ 2 \ Projecting above equation with (<fr vw \, 
we get the following expression for the correlation energy 

1 



1 



s 2 



\$ vw ) = AE 



corr 

VW ' 



(29) 



where we have used ($p,1,\& vw ) — 0. This equation can 
further be simplified, using Wick's theorem, as 



{H^S} + ~{H^S 2 }\$ VW ) 
E^(<S>Z\S+\s 2 \$ m 



(31) 



where E%™ is the difference between the exact energy of 
the closed-shell and two- valence states. It has the expres- 
sion 



patt 



(32) 



where, e v and e w are the Dirac-Fock energy of the va- 
lence orbitals \<f> v ) and \<j> w ) respectively. And AE%™ = 
AE^™ 1 — AEq ", is the difference of the correlation en- 
ergy of closed-shell and two- valence states. Diagrammat- 
ically, AE^ in Eq.(32) is equivalent to the closed dia- 



grams with free lines representing the valence states at 
the vertexes. Like in the second order MBPT, one can 
separate the AE^ diagrams to one- and two-body types. 
The one-body diagrams are similar to the ones in Fig. 
[T with the bottom interaction (dotted line) replaced by 
cluster amplitude (solid line). Similarly, the two-body 
diagrams are similar to those of Fig. [2jb-g) with the 
bottom interaction replaced by the cluster amplitude. 



C. CC equation from Bloch equation 

The CC equations discussed so far are derived from 
the eigenvalue equation of the Dirac-Coulomb Hamilto- 
nian. Another approach is based on the generalized Bloch 
equation given in Eq. ([5]). This is more transparent to 
implement and convenient to analyse the working equa- 
tions of CC with incomplete model space. In Eq. (|5j, 
the second term on the right hand side, renormalization 
term, is often defined as 



W = PVQP = (Vft)close. 



(33) 



Here, close indicates the operator connects states within 
the model space. Diagrammatically, the representation 
of the operator has no free lines in the closed-shell sector 
and only valence orbitals as free lines in the open-shell 
sector. Using Eq.(33l, we can write 



P,h }p = (yn-nw)P. 



(34) 
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Operating on the two-valence atomic reference state, 
|$i>.a;), and projecting with the doubly excited determi- 
nant |, we get 



V< 1 ( 1 + 5+ ^S 2 



1 + S+-S 2 )W 



I ^vw) ■ 



(35) 



From Wick's theorem further simplification follows af- 
ter contracting the operators. There are connected and 
disconnected terms, however, only the connected terms 



remain [T5] on both sides of Eq.(35). We get 



($%\{H S} - {SH }\* VW ) = 



Ve T l 1 



'-S' 



1 + 5+ -S 2 )W 



1$, 



5- 



(36) 



where the subscript conn refers to connected terms. To 
arrive at the equation we have used ($^|[T, H ]\& vw ) — 
0, as T being the closed-shell cluster operator, it does not 
operate in the valence space. To examine the equation in 
further detail, consider the terms on the right hand side. 
Expanding the exponential in the term Ve T 



= V 



{VT} 

m 



i F=n 

-{VTT} 



{VTTT} 



- — {VTTTT} = V, 



(37) 



where V is the dressed operator. Similarly, for the other 
terms 



e T SW) 
(e T W) 



{VS}, 
{SW}, 

o, 



(38) 

(39) 
(40) 



as no contraction can occur between T, the closed-shell 
cluster operator, and open-shell operator 5 to obtain con- 
nected term. The same is true of T and the effective 
interaction W. The reason is, T operates on the closed- 
shell sector, whereas W operates in the valence sector. 
Though it is not shown explicitly, there are similar re- 
lations for (l/2)5 2 as well. From the definition of the 
normal Hamiltonian Hh — V + Hq we can combine two 
of the terms as 



H S + VS = H N S. 



(41) 



Using Eqs. ( |37|40[ ) in Eq.@, we get the CC equation 
in the form 



($^|ff N + {tf N 5} 



l rr^ 

■-{SSH eS }\<f> vw ) = 0, 



(42) 



where H c g = Hq + W, is the effective Hamiltonian. The 
form of the effective Hamiltonian H e g is close, no free 



lines or only valence lines as free lines, therefore Eq.(42) 
can be written as 



1 



(V2l,\H N + {H N S} + -{H N S 2 }\$ VW ) 
H^ZlS+ls 2 ]^ 



(43) 



This is identical to Eq.(31), which is obtained from the 



eigen value equation of the Dirac-Coulomb Hamiltonian 
with the exponential ansatz. 



D. Diagonalization of H c ff 

In the single reference calculations the mapping from 
reference state to exact state is simple, and straight for- 
ward. The exact state, as given in Eq. Q, is the trans- 
formation of reference state \& vw ) with f2. It is not so 
simple with multi-reference model spaces. The model 
space then encompasses a set of determinantal states 
{l&VaWf))} € P, however, each state by themselves are 
not the reference states. The CC equation in Eq. ( 43 ) is 
then modified to 



<*5U I#n + {hJ} + ±{H N SS}\$ VaWll ) 



0,(44) 



7,<5 



where, the sum over S and 7 spans all the determinantal 
states within P. This is the working equation of multi- 
reference two- valence CCT with CMS. The last term re- 
quire careful consideration while implementing and as we 
mentioned earlier, folded diagrams arise from this term. 

The wave operator f2 is defined once the CC equations 
are solved, but the model functions are not yet defined. 
Next step of the calculation is then to evaluate the matrix 
elements of the effective Hamiltonian 



H eS (v, to; x, y) = (<f> vw \H + Vn\<f> a 



(45) 



The H c ff matrix is non symmetric as Q operates on the 
ket state and after diagonalization, one gets a biorthogo- 
nal set of eigen states l^) . These are the model functions 
of the multi-reference CC, the exact state is then 



l*i) = n|*?>, 

and the eigen value equation is 

The eigenstates in general are of the form 

i*°> = £<si*w> 

a/3 



(46) 
(47) 
(48) 



where c\^ are the coefficients of the linear combination 
or eigen vector elements of H c g . 



IV. INCOMPLETE MODEL SPACE 

Incomplete model space (IMS) consists of a restricted 
number of configurations from the CMS. Remaining con- 
figurations are part of the orthogonal space. Outcome 
of such a model space is, the clean separation of inter- 
nal and external cluster amplitudes is no longer true. 
Further more, the subsystem embedding condition is vi- 
olated. For example, cluster operators which are external 
in one-valence Hilbert space may no longer be so in the 
two-valence Hilbert space. The intermediate normaliza- 
tion Eq. (|6| is then, in general, not applicable 



(49) 



Following which, the H e g is not guaranteed to be opera- 
tional only within the model space, it may as well connect 
a state in P to a state in Q. Where as the obvious advan- 
tage of defining H c tf is to work within the model space 
and incorporate the effects of orthogonal space in an ef- 
fective way. Restoring the operational space of H c s to P 
requires a set of constraint equations [49] and a previous 
work reported the implementation of particle- hole sectors 
[52"] in relativistic CC calculations. However, all the good 
virtues of CMS, in the context of Fock-space CCT, are 
applicable when the model space is quasi-complete. For 
a lucid description of what constitutes a quasi-complete 
model space refer [J3J [5T] . 

Like in CMS, as a specific example consider the low- 
lying states of Yb. An ideal incomplete model space 
would consist of the configurations 6s 2 , 6s6p and 5d6s. 
Model space would then encompass all the levels im- 
portant to ongoing precision experiments: 6s 2 ( 1 5'o), 
6s6p ( 3 Pj), 5d6s ( 3 Dj) and 6s6p (^i). Obvious ad- 
vantage in such a selection of model space is isolation, 
as evident in Fig. [31 from the potential intruder states 
6p7p ( 3 -Pj), 6s7s ^Si) and 6s7s ( 1 5'o). Here, we can ap- 
ply subduction process to check if the model space con- 
sidered is quasi-complete and is shown in Fig. [6] Initial 



6s6p 




5d6s Two-valence 



5d y One-valence 



model space ( 6s, 6p and 5c?). Finally, removal of an- 
other electron gives the closed-shell model space. All 
the configurations obtained in the subduction are part of 
respective model spaces. This is a requirement of quasi- 
complete model space and necessary condition for sepa- 
ration of internal and external excitations. 



V. PROPERTIES CALCULATIONS 

A. Hyperfine structure constants 

The HFS constants of an atom are the parameters 
which measure further splitting of fine structure levels. It 
arises from the interaction of electromagnetic moments of 
the nucleus with the electromagnetic field of the atomic 
electrons |53j . The general form of the hyperfine interac- 
tion Hamiltonian is 



Hhts — 



= EE(- 1 Wi) T - 



(50) 



k . q 



where t£(r) and T* are irreducible tensor operators of 
rank k effective in the electron and nuclear spaces re- 
spectively. For the magnetic dipole interaction (k = 1), 
the explicit form of the tensor operators are 



tl(r) = 



-iv^[a-d(f)], 



cr* 



T 1 



(51) 



Here, Ci(f) is a rank one tensor operator in electron 
space and fi q is a component of /x, the nuclear magnetic 
moment operator. The HFS constants are the expecta- 
tion value of ifhfs and the magnetic dipole HFS constant 
is then 



(52) 



Closed-shell 



Where, \^f l ) is the exact wave function expressed in 
Eq.(17), using coupled-cluster theory. The denominator 
(^i|<F~) is the normalization factor and it is not be con- 
fused with intermediate normalization Eq. (|6| . The later 
defines the relation between the reference state and the 
exact state. And it does not determine the normalization 
of the exact state. 



B. HFS constant in one- valence sector 



FIG. 6: Incomplete model space of Yb two- valence calcula- 
tions. Arrows indicate the subduction to lower valence sectors 
and respective model spaces. 



stage is the two valence model space consisting of 6s 2 , 
6s6p and 5g?6s. Removal of one electron from each of 
the configurations leads to a configuration in one- valence 



Once the CC equations and cluster amplitudes are 
known, the atomic properties are calculated with the 
exact atomic states so obtained. It is expectation for 
dynamical variables and matrix element for transition 
amplitudes. From the CC wave function of one valence 
systems in Eq.(17), the expectation of -f/hfs is 



(*„|£Th&|tf„) = (^v\H M s + 2S^H hls + S^HM S S\^ v ), (53) 
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where, H his = e T ^ H h f s e T is the dressed operator. The 
factor of two in the second term on the right hand side 
accounts for H^S as S^-ffhfs = -HhfsS'- An expansion of 
i/hfs ideal for an order wise calculations is 



CO 

H hfs = H Ms e T + J2 -f (r f )" ^hf s e 3 



(54) 



n=l 



The normalization factor, denominator in Eq.(52), in 
terms of coupled-cluster wave function is 

(*„|*„) = ($,| (1 + S+) e Tt e T (1 + 5) (55) 

Note that dressed operator /f^fs and e T ^e T in the nor- 
malization factor are non terminating series. In a recent 
work, we demonstrated a scheme to include a class of 
diagrams to all order in T iteratively for properties cal- 
culations. With the method we calculated the magnetic 
dipole HFS constant of the singly ionized alkaline-Earth 
metals [18] . Based on the extensive calculations reported 
in ref . [18] , we conclude terms higher than quadratic in T 
contribute less than 0.1% to the HFS constants. So there 
are no compromises on important physics when Hm s , 
Eq.(57l, is truncated after the second-order in T. How- 



ever, there are enormous gains in computing resources 
and simplification of the procedure with the iterative 
scheme. Here, we shall not dwell on the iterative scheme, 
interested readers may refer to ref. |18j for more details. 



C. HFS constant in two- valence systems 



From the CC wave functions of two-valence systems 



defined in Eq. ( 26 1 , we get 



S + -S 2 ) Ih 



Ill's 



S + -S 2 



//,„., ( S+ -S 2 ) hi.,,} 



(56) 



Where to shorten the notations we have replaced the va- 
lence indexes in the two- valence states v a wp (vsw 7 ) with 
j (k). This is the CC expression to calculate the HFS 
constants of two-valence electron atoms. The operator 
H^fs is, as defined earlier, the dressed HFS interaction 
Hamiltonian. As discussed in the one- valence case, com- 
prehensive inclusion of all order of T is beyond the scope 
of current theories and computational resources. Hence, 
for the two-valence sector we consider up to quadratic 
terms of T in Hm s , approximately 

H Ms « H hls + H his T + T^Hms + T^H his T. (57) 

To compute -ffhfs of the two-valence sector, we bor- 
row the concept of effective one- and two-body operators 



from our previous work ref [18] . Diagrammatic repre- 
sentation of the effective operators are as shown in Fig. 
[7] It is important to note that the two-body effective 
operator, shown in Fig. [7|e), arises from the last term 
in Eq. (57). And, since it has two-orders of T, the ac- 



tual hyperfine diagrams obtained after contraction with 
S may have negligible contributions. For this reason we 
shall not elaborate on the HFS diagrams arising from the 
dressed two-body effective properties operator. However, 
we do incorporate these diagrams in the calculations and 
mention the contributions in the results. The diagrams of 



(a) 



(b) 



(<0 



(d) 



(e) 



FIG. 7: Representation of effective one- and two-body dressed 
properties operators. 



the HFS constant in the two-valence sector are grouped 
into different categories. First based on the cluster oper- 
ators in the expression and later, in terms of the number 
of core, valence and virtual orbitals. Specific terms and 
groups of diagrams are discussed in this work. 



1. Effective one-body operator 

There are four diagrams from -ffhfs which has non-zero 
contribution. These are the two- valence diagrams from 

T 1 

T2 with the bare hyperfine interaction h^{ s inserted to 
all the possible orbital lines. Contributions from these 
diagrams is expected to be very small as S are not a part 
of the diagrams. 




(T 




(d) 



(f) 



FIG. 8: Hyperfine diagrams contributing to the terms, 
Si^ H his (diagram (a)), S (2)t H his (diagram (6)), 



^-ffhfs^ 23 (d 

(diagram (/)). 



S. 



lagrams 



form (c) to (e)), and S (2) H hfs S {2) 



Diagrammatically, Hy 1 { s S < - 1 " > and S^ 1 ^ each have 
one diagram and these arise when is contracted with 
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the one-body effective operators: with diagram in 

Fig. 7[d), and with diagram in Fig. 7fc). The 



diagram arising from ffhfs is as shown in Fig. [bj^a J 
Time reversed version of the same diagram correspond to 
HhfsS^-, however, this is not shown in figure. The con- 
tributions from Hhf s S^ and .ffhfs are large as these 
are only first order in S. Further more, ffhfs is one-body 
interaction and hence, one-body effective interaction are 
dominant. 

The terms HmsS^ and S^ 2 ) -ffhfs each have one dia- 

(2) 

gram and These arise from the contraction of Sv; with 
the one-body effective operator of ffhfs shown in diagram 
Fig. [7'a). The diagram from S^^ffhfa is shown in Fig. 



|8|Jb) . Like in the previous case, the time reversed diagram 
arise from ffhf s <S^ 2 ' and is not shown in the figure. One 
can expect these terms to constitute the leading order 
as these are the lowest order terms with S% . Rationale 
for such an anticipation is, in general, the magnitudes of 



("2") 1 . ! < 

are larger than S^' and T operators 



(i) 



(a) 



(b) 



(c) 



(<1) 



FIG. 9: Diagrams arising in contraction of with Srp. 



2. S r H h{s S 



The leading term in Eq. ( 56 ) which is quadratic in S 



S^ffhfs^+C.C 



(58) 



where, c.c. represents complex conjugation. We now 
discuss the diagrams arising from each of these terms. 

There are sixteen diagram arising from S^^ffhfsSW and 
topologically, these are the effective one-body diagrams 

Fig.^a-b) sandwiched between and S^. To exam- 
ine the diagrams in more detail, all the diagrams (four in 

all) from the contraction s!^~Srp are shown in in Fig. [9] 
To each of the diagrams in Fig. [9] the effective one-body 
operator can be inserted in four ways. As an example, 
consider the diagram in Fig. [9jb) , all the four diagrams 
after inserting the effective one-body operator are shown 
in Fig. [10] 

There are five diagrams from SW T ff h fsS" (2) . These 
arise from the contraction of with through one- 
body operators in Fig. (7ja) and (c). The diagrams from 



(d) 



FIG. 10: Hyperfine diagrams obtained after inserting the Hy- 
perfine interaction operator in diagram (b) of Fig.|9|. 



£>2 "Ms 



are as shown in Fig. |8jc-e). However, the 
diagrams from are not shown. Identical number of 

diagrams arise fromS , ( 2 ) T ff h fsS' (1) . The effective diagram 
in this case are Fig. [7ja) and (d). 

Finally, only one diagram arises from the last term, 
S^H h{s S^. This diagram is shown Fig. j^f). Only 
Fig. [7]^a) is the allowed effective one-body operator which 
contribute to this term. 



D. Electric dipole transition amplitudes 

Electric dipole is the most dominant electromagnetic 
multipole in the radiative transition of atoms. In ma- 
jority of the cases, depending on the decay channels, it 
defines the life time of an excited state. Theoretically, 
the relevant quantity is the reduced matrix element of 
the dipole operator D between the initial and final states 
and |\&/), respectively. The two states are opposite 
in parity as D is an odd parity operator. The expression 
of the reduced matrix element is 



(59) 



Once the reduced matrix elements are evaluated, the ac- 
tual matrix elements of the specific states are calculated 
from the Wigner-Eckert theorem. Here, we need to make 
a finer distinction of the wave operator in the valence 
universal or Fock-space CCT. As ff DC commutes with 
parity, so does the wave operator Q, and we can consider 
f2 as 



(60) 



Where f2 + and operates on the even and odd parity 
reference states within the model space. The separation 
of these two components follows naturally from the parity 
selection rules imposed on the cluster amplitudes. How- 
ever, this ought to be handled with care as complications 
arise in the calculations of peturbed cluster amplitudes. 
Where there is a second perturbation, besides the resid- 
ual Coulomb interaction, which is parity odd. We may 
rewrite Eq. ( 59 ) more precisely as 



D 



if 



jWfWfWW 



(61) 
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In the one-valence sector, the reduced matrix element of 
D is 

(y w \\D\\V v ) = (9 V \\D + S^D + DS + S^DS\\9 V ), (62) 

where, \9 V ) and \^ w ) are the initial and final states in 
terms of the valence states. Though the expressions are 



similar to Eq.(53), there arc two important differences. 



Unlike the HFS energy splitting expression S'D ^ DS, 
this is because the initial and final states are different. 
Same set of properties diagrams in HFS calculations, af- 
ter modifications to account for the two key differences, 
are then adopted to compute reduced D matrix elements 
with CC wave functions. 

After a similar modification, like in Eq. ( 56 1 , for the 
two- valence systems 



;<!>, D + D[ s+ l -s 2 



S+^S 2 j D+ [S+-S 2 



D[S+-S 2 \ \<S> k 



(63) 



Where the notations and terms are the same as in the 
two-valence HFS case. However, the two key differences 
mentioned earlier still hold true. Like in HFS 



D + DT + T 1 D + T^DT. 



(64) 



We then proceed like in HFS and calculate the effective 
diagrams, both one- and two-body. These are then con- 
tracted with the cluster operators and we evaluate the 
transition matrix. 



VI. NUMERICAL DETAILS 

Coupled-cluster theory, though powerful is computa- 
tionally intensive and implementation is non trivial. The 
large number of unknowns and equations demand special 
attention to all aspects of computations. Right from the 
initial stage of identifying and calculating the cluster di- 
agrams, to the final stages of solving the CC equations 
and computing properties from the CC wave functions. 
Here, we give concise description of what we consider ab- 
solutely essential, theoretical and computational aspects 
for atomic coupled-cluster calculations. This choice is en- 
tirely based on our experience of developing and imple- 
menting CCT spanning closed-shell, one- and two- valence 
systems. Besides the CC wave-function calculations, we 
have also proposed, developed and implemented methods 
to compute properties from CC wave functions. The se- 
lected issues addressed are provided with the anticipation 
that interested researchers shall find these details valu- 
able. And facilitate adopting CCT for atomic many-body 
computations with minimal effort. 



A. Orbitals and basis functions 

Results presented in this paper are based on the Dirac- 
Coulomb Hamiltonian H DC given in Eq. ([lj. It incor- 
porates relativity at the single particle level accurately. 
And, as the name indicates, the Coulomb interactions 
between the electrons. For the nuclear potential Vm{t) : 
we consider the finite size Fermi density distribution 



:(r) 



Po 



1 + e( r_c )/ a ' 



(65) 



here, a = £4 In 3. The parameter c is the half-charge 
radius, that is p n uc(c) = po/2 and t is the skin thickness. 
At the single particle level, the spin orbitals are of the 
form 

where P nK (r) and Q nK (r) are the large and small compo- 
nent radial wave functions, k is the relativistic total an- 
gular momentum quantum number and X«;m( r / r ) are the 
spin or spherical harmonics. One representation of the 
radial components is to define these as linear combination 
of Gaussian like functions and are referred as Gaussian 
type orbitals (GTOs). Then, the large and small compo- 
nents EH [55] are 



P„ K (r) 



(67) 



The index p varies over the number of the basis functions. 
For large component we choose 

g L Kp {T)^C L Kl r^e'^\ (68) 

here n K is an integer. Similarly, the small component is 
derived from the large components using kinetic balance 
condition |56j . The exponents in the above expression 
follow the general relation 



(69) 



The parameters ao and (3 are optimized for each of the 
ions or neutral atoms to provide a good description of the 
properties. In our case the optimization criteria are to 
reproduce the numerical result of the self consisten field 
(SCF) energy and orbital energies. 

From Eq. ( 66 1 the reduced matrix element of the mag- 
netic hyperfine operator between two spin orbitals , v' 
and v, is 



(l/H^Hf) = — (k v + K v ')( — K v ' \\C\\n v ) X 

dv 

A detailed derivation is given in Ref. [57] , 



K ,p nvK jy.o) 
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For the alkaline Earth metal atoms Sr and Ba as well 
as Yb, we use V N ~ 2 orbitals for the calculations. This is 
equivalent to calculating the spin orbitals from the sin- 
gle particle eigenvalue equations of the doubly ionized 
atoms, namely Sr 2+ , Ba 2+ and Yb 2+ . The single par- 
ticle basis sets then have few bound states and rest are 
continuum. The basis set is optimizeed such that: single 
particle energies of the core and valence orbitals are in 
good agreement with the numerical results. For this we 
use GRASP92 [58] to generate the numerical results. 



B. Orbital subsets 

Orbitals, the single electron wave functions, in closed- 
shell systems are separated into two distinct subsets, core 
(occupied) and virtual (unoccupied). Former are shells 
which are completely filled in the ground state determi- 
nantal state and later are empty. Distinction is not so 
straight forward in open-shell systems. The classification 
of the orbitals for Yb atom in our current calculations 
is shown in Fig. 
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The valence orbitals are partially 
filled in the model functions and are like core orbitals, 
electrons can be excited from the valence shells. This is 



Virtual 




Is 



IF 



Valence 



Core 



FIG. 11: Classification of orbitals into core, valence and vir- 
tual subsets. Few orbitals are members of more than one 
subset. 



particularly true when considering the valence as parti- 
cles. Consequently, as discussed in the next subsection, 
the closed-shell diagrams can be modified to the open 
shell ones. On the other hand, valence shells can also ac- 
commodate excitations from the core shells. A property 
typical of shells in the virtual space. Hence, in the cluster 
amplitudes the excited states incorporate the valence or- 
bitals as well. The dual character of the valence orbitals 
can be adapted for faster diagram evaluations with ap- 
propriate rearrangement of the summation sequence. For 
example, there is enormous computational advantage in 
considering the free orbital lines, in the cluster diagrams, 
as the outer most summation. The example given and 



many other features we have developed are more compu- 
tational in nature and less of physics. We shall elaborate 
on these matter in future publications devoted to the 
computational aspects of our work. 



C. CC equations 

In Fock-space CCT, as mentioned earlier, the cluster 
operators are generated sector wise in sequence. First, 
the close d-shell cluster amplitudes are computed from 
the Eq. ( 20j21 ). Where the dressed operator J$n in the 
closed-shell CC equations, like in Eq. (371, is 



{#nT} 



i F=n 



i p^n 

-{tf N TTT} 



--{H N TTTT}. 



(71) 



The closed shell CC equations as evident from the expres- 
sion of Hn arc nonlinear equation. In CCSD approxima- 
tion, we have second and fourth order nonlinearities in 
T 2 and Ti, respectively. However, the working equations 
can be written in linear form as 



A u (T)Ti+Ai2(T)r 2 - B u 
A 21 {T)T 1 + A 22 (T)T 2 = B 2 . 



(72) 
(73) 



Since the original equations are nonlinear equations, the 
coefficients Aij (T) are functions of cluster amplitudes T. 
On the right side B are the matrix elements of H^. The 
equations are then solved iteratively till convergence. 



(a) 



0>) 



(c) 



FIG. 12: Conversion from closed-shell cluster operator T 
diagrams to open shell operators S. 



To set up the equations, we evaluate the terms based 
on diagrammatic analysis. There are several diagrams 
and for example, those arising from linear terms are given 
in ref. [TB]- The total number of equations scale as 
N^Nq, where N v and N a are number of virtual and occu- 
pied orbitals, respectively. For the calculations discussed 
in this paper N v 130 or more and N Q w 20. The 
coefficient matrix A is non-symmetric and dense, so the 
number of matrix elements scales as ~ N^N^, which is 
~ 4.6 x 10 13 for typical examples in the present compu- 
tations. It is an extremely large matrix and impractical 
to store. In addition, the elements are functions of T 
and not static. For these reasons, the elements of A are 
calculated on the fly, as and when needed. Another com- 
plication is, the operations required to generate each ele- 
ment of A in the equations scale as N^N 2 . All together, 
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combining the number of matrix elements and number 
of operations, number of binary arithmetic operation in 
each iteration is 0(N®N®). Which is indeed a very large 
number for high Z atoms. 

Diagrammatically, to generate the closed-shell equa- 
tions, we identify all the diagrams in the closed-shell CC 
equations and evaluate the angular integrations based on 
angular momentum diagrams |16) . An example diagram 
is shown in Fig. [l2[a), it is the double contraction of 
V with T2 and contributes to the T2 equation. To set 
up the one-valence and two-valence cluster equations we 
avoid diagrammatic evaluation. Instead, the closed-shell 
diagrams are topologically transformed into one-valence 
diagrams. As shown in Fig. 12 'b), one of the core or- 



bital line is rotated and transformed into a valence line. 



Diagrams so obtained are very different from Fig. 12 a) 
in terms of possible contractions. However, the results 
from the angular integration remain unchanged. A sim- 
ilar procedure is adopted for the two-valence equations 
as well. With this, a careful analysis and evaluation of 
closed-shell CC equations is the only requirement to set 
up the one- and two-valence CC equations. The coupled 
nonlinear and linear equations are solved iteratively. We 
employ direct inversion in the iterated subspace (DIIS) 
[59] for convergence acceleration. 



VII. RESULTS 

A. Excitation Energies 

Excitation energies of the one-valence ions Sr + and 
Ba + were reported in our previous paper |18j . Here, we 
report the results of the calculations in the two-valence 
sector, the excitation energy of an state nln'l' ( 2S+1 ^Lj 



, from Eq. (47 1, is 



AE. 



nln'l' ( 2S +i-)Lj 



E. 



nln'l' (2S+l) L/ 



E„ 



(74) 



Where, E n 



'-So 



and E nln , v (2 s+i)i / are the ground and 



(47). Vaeck, 



excited state eigenvalues of H e g in Eq, 
Godfroid and Hansen [BOJ had calculated t 



le excitation 

energies and investigated in detail the configuration mix- 
ing of atomic Sr with multiconfiguration Hartree-Fock 
theory with special attention on the singlet states 1 L. 
In particular the states with configurations of the form 
5snp 1 P°, 5snd X D and bsnf 1 F°, including the Ry de- 
berg states. Improved experimental data and prospects 
of cooling and trapping had spurred further theoretical 
studies on properties of Sr. Important recent theoretical 
work are by Porsev and collaborators [61] . and Savukov 
and Johnson [62]. Previous works of Eliav, Kaldor and 
Ishikawa [TO] IM] reported the excited energy calcula- 
tions of atomic Ba and Yb using the Fock-space based 
coupled-cluster theory. The other widely used atomic 
many-body method employed is CI-MBPT [65], based 
on Dzuba and Ginges [30] calculated the excitation en- 
ergies of Ba. Same method was used by Porsev and col- 



laborators to calculate the excitations energies and HFS 
constants of Yb [67] . 

One reason for choosing the three atoms in our cal- 
culations is the significant difference in the sequences of 
ns(n — l)d 3 Dj, ns(n — l)d 1 Dn^nsnp 3 Pj and nsnp 1 Pi 
levels. As evident from Table. |lj in Sr the 5s4d 2S+1 Dj 
levels lies between 5s5p 3 Pj and 5s5p 1 P\. Whereas the 
6s5d 2S+1 Dj levels are below the 6s6p 2S + 1 p J levels 
in Ba. The difference in the level structure can be at- 
tributed to the presence of an additional diffuse shell Ad. 
The sequence gets more complicated in Yb, 6s6p 3 Pj 
levels are below 6s5d 3 Dj, however, the 6s6p 1 P\ lies 
between 6s5d 3 Z?2 and 6s5d 3 D%. It is to be noted that, 
the difference between Ba and Yb configurations is the 
presence of 4/ in the Yb core. And, is the cause for the 
change in the level sequence. 

The excitation energies obtained from our calculations 
are reasonably close to the other theoretical results for 
Sr and Ba. However, there is a lack of clear trend in 
the differences. For the excitation energies of Sr, our re- 
sults are consistently better than the MCHF results [BOJ . 
And, our result of 5s5p 3 P2 is closest to the experimental 
value. One observation is, although Porsev and collabo- 
rators [61] , and Savukov and Johnson |62j used the same 
method CI-MBPT, the results from the former are con- 
sistently better than the later. Possible reason could be 
the single particle basis set. The former used a combi- 
nation of V N , V^ -1 and V N ~ 2 orbitals for the core and 
valence, and virtuals are generated through a recurrent 
procedure. In the later work, the orbital set are B-splines. 
As described earlier, we use numerical Gaussian type or- 
bitals for our calculations. 

Comparison of excitation energies of Ba presents an 
interesting case. Some of our results are better than the 
previous CC results of Eliav and collaborators [M] . The 
results provides a numerical validation of our approach. 
On the other hand, the results of Dzuba and collaborators 
[65] . and Safronova and collaborators [66] uses similar 
basis sets but different many-body methods. The former 
used CI-MBPT, whereas the later used the recently de- 
veloped CI plus all order method [66]. The results from 
the later are consistently better than the former. 

The Yb excitation energy calculations presents a seri- 
ous challenge. Earlier CC calculations of Yb excitation 
energies [19] could not reproduced the experimental se- 
quence and we also encounter the same issue. In partic- 
ular, the 6s6p 1 P\ is above the 6s5<i 3 Dj levels, whereas 
experimentally it lies between 6s5d 3 £>2 and 6s5d 3 D%. 
Sequence in our results is similar and individual values 
are consistently higher than the previous CC calcula- 
tions. The sequence, however, is correctly reproduced 
in another calculation with the CI-MBPT method [67] . 
where the basis set used is combination of Dirac-Fock or- 
bitals for the core and valence, and virtuals are generated 
through recurrent procedure. The comparison of the dif- 
ferent results indicate a wide variation in the many-body 
methods and single particle states. In fact, none of the 
works listed and referred have a common many-body the- 
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TABLE I: Two-electron removal energy and the excitation energies calculated using relativistic coupled-cluster theory. All 
values are in atomic units. 



OlclLc 




W III ItioLllL 


Other work 


.cjA-jj icoliil x\ei IUOI . 






EE 




EE 




EE 






Atnmir 87c !r- rKrl t U 2 










5s 2 


— n fii 93Q 


0.0 




0.0 




0.0 


5s5p 3 P 


-0.55169 


0.06771 




0.06566 a 




0.06524 


5s5p 3 Pi 


-0.55170 


0.06768 


0.06651 a , 0.06871 b 


0.06609 


K o E 3 p 

DSop tr% 


— U.OOZUo 


U.UD ( oo 




0.06833 a 




0.06788 




— U.OOOOl 


U.Uoooo 




0.08230 a 




0.08274 


KcAfJ 3 Do 


U. iJO' L ± 1 o 


084fi1 

U.UOIUl 




0.08260 a 




0.08301 


5<i4r7 3 Do 








0.08312 a 




0.08347 


5 s4<7 1 Do 

'. 7 f ' 1 ( 4, ' - _ 


—0 52594 


09345 


0.09210 a , 0.11477 c 


0.09181 


5 s4rf 1 Pi 


—0 51 283 


1 0656 


n noQ^i a 
U.Uyool , 


U.lUUlo , 


0.10730 c 0.09887 






Atomic 137 Ba; [Xe]6s 2 










6s 2 X S 


-0.56439 


0.0 




0.0 




0.0 


6s5d 3 Di 


-0.52303 


0.04136 


0.04211 d 


0.04106 c , 


0.04119 f 0.04116 


6s5d 3 D 2 


-0.52170 


0.04269 


0.04296 d 


0.04193 c , 


0.04200 f 0.04199 


6s5d 3 L> 3 


-0.51960 


0.04479 


0.04473 d 


0.04375 c , 


0.04366 f 0.04375 


6s5d x Do 


—0 51030 


05409 

\ f • \f J IV/-/ 


0.05395 d 


0.05197 c , 


0.05298 f 0.05192 


6<s6« 3 Pn 


—0 50667 


05772 


0.05697 d 


0.05575 c , 


0.05591 f 0.05589 


6<s6« 3 Pi 


—0 50540 


05899 


0.05869 d 


0. 05742°, 


0.05758 f 0.05758 


6s6» 3 Po 


—0 50311 


06128 


0.06284 d 


0.06147 c , 


0.06159 f 0.06158 


6s6n x Pi 


—0.47291 


09148 


0.08409 d 


0.08256 e , 


0.08125 f 0.08229 






Atomic 173 Yb; [Xe]4/ 14 6s 2 










6s 2 


-0.68083 


0.0 




0.0 




0.0 


6s6p 3 P 


-0.59944 


0.08140 


0.07909 g , 


0.07874 h , 


0.07877 


0.07877 


6s6p 3 Pi 


-0.59645 


0.08439 


0.08242 g 


0.08200 h , 


0.08200 


0.08198 


6s6p 3 P 2 


-0.58914 


0.09170 


0.09038 g 


0.08999 11 , 


0.09002 


0.08981 


6s5d 3 Di 


-0.56110 


0.11973 


0.11362 g , 


0.11425 h , 


0.11158 


0.11158 


6s5d 3 D 2 


-0.55975 


0.12109 


0.1 1473 s 


0.11136 h , 


0.11274 


0.11278 


6s5d 3 D 3 


-0.55602 


0.12481 


0.11699 g , 


0.11503 h , 


0.11517 


0.11514 


6s6p x Pi 


-0.55301 


0.12782 


0.12426 g , 


0.11253 h , 


0.11669 


0.11422 


6s5d 1 D 2 


-0.54667 


0.13416 


0.13025 s 


0.12595 h , 


0.12672 


0.12611 



a Reference [51]. b Reference [52]. 

c Reference [50]. d Reference [53]. 

c Reference[66 . f Reference 30 . 

s Reference [19]. h Reference [67] • 
1 Reference [68] . 



ory and basis sets. Perhaps, this is an indication of the 
issues which require consistent efforts to resolve the dif- 
ficulties of precision calculations of two-valence systems. 

B. Hyperfine structure constants 

Hyperfine constants are appropriate atomic properties 
to inspect the accuracy of atomic wave functions in the 
small radial distances-within and close to the nucleus. 
For the calculations we use the nuclear moments given 
in the compilation of Stone [69) and the values are given 



Table, [n] 



TABLE II: The nuclear spin I, the magnetic moment fi (in 
nuclear magneton), and the electric quadrupole moment Q 
(in barn), used in the paper. 



Ion 


I 




Q 


87 Sr 


9/2 


-1.0936(13) 


+0.305 


137 Ba 


3/2 


+0.93737(2) 


+0.246(2) 


i73 Yb 


5/2 


-0.648(3) 


+2.80(4) 
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Hyperfine constant A of Yb + 



3. Hyperfine structure constant B of d states 



The magnetic dipole HFS constant of 173 Yb + and the 
electric quadrupole HFS constants of 87 Sr + , 137 Ba + and 
173 Yb + from our calculations are given in Table 



III 



For comparison, the results from other theoretical stud- 
ies and experimental data are also given. Contributions 
from the specific terms in the CC properties expression 
of HFS are listed in Table. |IV| Previous theoretical study 
by Martensson [71] on the magnetic dipole HFS con- 
stant of 173 Yb + is based on the CCT and basis set is 
obtained from discrete spectrum method [70]. In this 
work, the HFS constant a of 6S1/2, 6pi/2 an d 6p 3 / 2 are 
calculated. Our result of 6si/ 2 is slightly higher than 
both the theoretical and experimental values. However, 
for 6pi/2 our result is lower than the previous theoret- 
ical result of Martensson [7T] and closer to the experi- 
mental data [7T]. Similarly, our result of 6^3/2 is lower 
than the value of Martensson [71]. Although, the many- 
body methods employed in the two calculations are the 
same, HFS constants of p in our results are lower. Com- 
ing to the 5d states, the previous calculations of Itano 
is based on multiconfiguration Dirac-Fock (MCDF) [58] , 
The 5d 3 / 2 results are close to our value, however, for the 
5tfg/2 state our results are much larger. At this stage it 
is difficult to pinpoint the reason for the large discrep- 
ancy between the two results. One observation from the 
component wise contribution in Table. |IV| is the large 
cancellation between the Dirac-Fock (leading order) and 
the next to leading order (S'Hhls + c.c). Similar pat- 
tern is observed in the HFS constant of nds/2 state of 
all the alkaline-Earth metal ions reported in our previ- 
ous work [18] . A comparison shows the cancellation is 
larger in Yb + . Another notable difference in Yb + is, the 
(S^HSi + c.c.) is large and almost cancels with S\HS2- 



The HFS constant of Sr + has been studied in sev- 
eral theoretical works. Our results of B are systemat- 
ically lower than the other theoretical values, which is 
evident from Table. |III| The previous calculations of 
Martensson [72] used the same many-body method like 
ours, relativistic coupled-cluster, but a different type of 
single particle basis set. The calculations of Sahoo [76], 
in terms of theoretical methods, is closest to ours. They 
have used the relativistic coupled-cluster and gaussian 
type basis set like we have done. However, Sahoo calcu- 
lated only for the 4(1^/2 state and his result is closest to 
the experimental data. Ours on the other hand is ~ 2.3% 
lower than his result. Similarly, for the same reasons for 
the HFS constants B of Ba + , the previous calculations 
of Sahoo [78] is closest to ours. However, our result of 
5d 3 /2 is closer to the experimental data. For 5d 5 / 2 state 
our value is ~ 0.6% lower than the theoretical value of 
Sahoo and ~ 5.3% higher than the experimental value. 
Considering that for the Sr + and Ba + calculations, the 
many-body method and type of single particle basis we 
have used are the same as in ref. [75] and [75], respec- 
tively. The difference in the results could be on account 
of minor differences in the exponents used in the basis 
set generation or the truncation of the coupled-cluster 
properties expression. There are striking changes, when 
compared with the P3/2, in the component wise contri- 
bution. Dirac-Fock contribution in both, Sr + and Ba + , 
are approximately three times larger than (S^H^fs + c.c. 
In addition, the contribution from the -f/hfs — DF, which 
essentially arises from the closed-shell part, is relatively 
large. This could be due to the diffuse electron density 
of the d orbitals and hence stronger interaction with the 
core electrons. 



2. Hyperfine structure constant b of p 3 / 2 states 



Leading terms, listed in Table. |IV| are the Dirac- 
Fock and (S^Hms + c.c). The later subsumes the core- 
polarization effects. For all the ions, Sr + , Ba + and YB + , 
the contributions from these two terms are almost equal. 
This is a significant deviation from the observed pattern 
in the magnetic dipole H FS constant [TSJ, which notice- 
able for Yb + in Table. IV Among all the theoretical 
calculations our results for Sr + is in better agreement 
with the experimental data. For Ba + , no previous theo- 
retical works and experimental data are available. Ours 
is the first study on the electric quadrupole HFS constant 
of the 6^3/2 state. Our result of Yb + is higher than the 
previous theoretical results of Martensson-Pendrill and 
collaborators [7T] as well as the experimental results of 
Berends and Maleki [531. 



Unlike the other two ions, Yb + has not been studied 
in fine detail. The previous theoretical work of Itano [75] 
is based on the MCDF method. And there are no ex- 
perimental data available for the 173 Yb + isotope. Our 
results are lower but close to the values from Itano [72"] . 
A closer inspection of the results from Itano's calcula- 
tions for the other ions (Sr + and Ba + ) reveals that, his 
results are consistently higher than the other theoretical 
and experimental data. One possible reason could be the 
contracted nature of the virtual orbitals, referred as cor- 
relation orbitals, in MCDF calculations. Hence we can 
expect a similar trend in Yb + as well and it is possible 
that our results are closer to the actual values. Com- 
pared to Sr + and Ba + , there is one remarkable change 
in the component wise contribution. There is large con- 
tribution from (S'Hhfs + c.c), which implies that there 
core-polarization effect is very important. It is on par 
with the Dirac-Fock term. 
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TABLE III: The magnetic dipole HFS constant for 173 
173 Yb+ ions. All the values are in MHz. 


Yb + and the electric quadrupole HFS constant for 87 Sr + , 137 Ba + , and 


Ion 


state 


This work 


Other works 


Experiment 


173 Yb+ 


6S1/2 
6pi/2 
6^3/2 
54/2 

54/2 


Magnetic dipole HFS constant A 
-3529.660 
-612.362 
-88.973 
-104.479 
22.078 


-3507 a 
-638 a 
-107 a 
-110. 31 b 
3.47 b 


-3497.5(6) a , -3508(9) c 
-518.2(4) a , -600 c 

- 


87 Sr+ 


5P3/2 

44/2 
44/2 


Electric quadrupole HFS constant B 
84.806 
33.961 
48.055 


82.655 d , 83.662 c 
35.075 d , 36.051 6 , 39.60 b 
48.800 d , 51.698 6 , 56.451 b 49.166 f 


88.5(5.4) s 
49.11(6) h 


137 Ba+ 


6P3/2 

54/2 

54/2 


98.954 
45.765 
62.685 


92.275 1 
51.32 b , 47.3 k , 46. 82 1 
68.16 b , 63.2 k , 62.27 1 


92.5(0.2) j 
44.541(17) k 
59.533(43) k , 60.7(10)', 62.5(40) m 


173 Yb + 


6P3/2 
54/2 
54/2 


1839.779 
902.301 
1165.046 


1780 a 
951. 4 b 
1190.4 b 


1460(50) n 


a Reference [71] . 
c Reference [75] . 
c Reference [75]. 
s Reference [77]. 
1 Reference |78|. 
k Reference [80]. 
m Reference 1521. 


b Reference [72]. 
d Reference [73] . 
f Reference |76|. 

Reference [38]. 
j Reference [79], 
1 Reference [51]. 
n Reference [551. 







4- Two-valence 

There are few theoretical and experimental work on 
the HFS constants of the neutral alkaline-Earth metal 
atoms and Yb. However, the importantance of such in- 
vestigations are likely grow in the near future as these, 
in particular Sr and Yb, are candidates of precision ex- 
periments and have been cooled to quantum degeneracy. 
In our work we make an effort to understand the sys- 
tematics to initiate a deeper analysis on the role of the 
electron correlation effects to properties like hypcrfine. 
The previous theoretical calculations of 137 Ba [53] and 
173 Yb [B7J are based on the CI-MBPT method and basis 
with different central potentials. For this reason it is non 
trivial to comment on the role of the correlation effects in 
a precise manner through a comparative study. The re- 
sults from our calculations, along with the leading order 
contributions, are listed in Table. |Vj From the table it is 
clear that, in most of the cases our theoretical results are 
not in very good agreement with the experimental data. 
Origin of the discrepancy could be the nature of the sin- 
gle particle basis we have used, the V N ~ 2 orbitals. On 
account of the doubly ionized charged state of the core, 
the orbitals are highly contracted and interacts rather 



strongly with the nucleus. Such orbitals are suitable for 
properties calculations of singly ionized states but not 
ideal for the neutral atoms. 

A very important aspect of our present work is the ob- 
served trend in the contributions from various terms. It 
is evident from Table. [Vjthe DF contribution is signifi- 
cantly dominant, it is far larger than the next to leading 
order contribution from what we refer as the one-body 
terms. The details of the one-body terms are discussed 
in Section. |V C 1| To quantify the relative contributions, 
define 

one-body terms 



Essentially the ratio between the leading and next to 
leading order contributions. The dominance of DF is par- 
ticularly true in the case of the ns(n — l)d 3 Dj states, 
among these states highest g is « 0.12 ( 5s4d 3 £>2 state 
of 87 Sr). For the nsnp 3 Pj states, the contribution from 
the one-body terms is small but not negligible. Largest 
and smallest value of g for these states are rj 0.3 for the 
6s6p 3 Pi state of Ba and « 0.23 for the 5s5p 3 P2 state 
of Sr, respectively. Other states have g close to 0.25. 

For the singlet states nsnp 1 P 1 and ns(n — l)d 1 -D 2 , 
the deviations from the experimental data are very large. 
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TABLE IV: Magnetic dipole and electric quadrupole HFS constants contributions from different terms. 



Ion 


state 






Coupled-cluster terms 








DF 


i?hf s -DF 


-t/hfs SlHhlaSl 
+C.C. +C.C. 


SIHmsSi Sl&hfsS-z Other terms Norm 



HFS constant A 



3 Yb+ 



7 Sr+ 



7 Ba 4 



3 Yb+ 



6S1/2 


-2582.096 


130.175 


-998.855 


-31.463 


-48.566 


-60.815 


14.667 


1.013 


6P1/2 


/ins ^ofi 


lo.oyo 




— / .UD4 


— iu.y / z 


— 0.101 


A QVG 
— 4.0 / O 


i on 


6P3/2 




1 voc 
l. /yo 


QO A 1 Q 


1 81 7 

— l.ol i 


1 OQyt 
— 1.ZD4 


*7 A A £ 
— I .440 


— U.oyl 


i on 
l.Ull 




—75.876 


0.192 


—19.903 


—0.870 


—1.289 


—8.127 


—0.122 


1.015 


5^5/2 


-28.899 


-0.536 


52.927 


4.081 


-0.424 


-4.903 


0.082 


1.011 








HFS constant B/Q 












6P3/2 


166.993 


-2.509 


109.511 


4.024 


2.565 


7.769 


-2.302 


1.001 


5^3/2 


80.939 


7.879 


26.888 


-.573 


1.073 


-3.540 


-0.791 


1.005 


5d 5 / 2 


110.863 


15.702 


34.180 


-.848 


1.377 


-1.911 


-1.081 


1.005 


6P3/2 


229.303 


-5.962 


170.985 


7.425 


5.364 


-1.153 


-3.097 


1.002 


54/2 


135.098 


12.635 


46.886 


.211 


1.305 


-7.849 


-1.213 


1.006 


5d 5 / 2 


172.975 


25.621 


63.015 


.333 


1.627 


-5.898 


-1.508 


1.005 


6P3/2 


372.894 


-16.855 


274.025 


11.547 


9.994 


-0.212 


6.441 


1.001 


54/2 


199.032 


11.143 


112.382 


1.292 


3.058 


-4.128 


0.975 


1.005 


54/2 


234.438 


21.220 


152.106 


2.217 


3.080 


4.774 


-0.041 


1.004 



A similar trend was also observed in the case of the ex- 
citation energy of these states as well. 



C. El transition amplitude 

We calculate the reduced matrix element of the dipole 
operator D from the expression given in Eq. (59). Once 



again, like in the HFS constants, we calculate the re- 
duced matrix elements of the Sr + , Ba + and Yb + ions 
from the intermediate one-valence wave functions. In 
present work, we do not attempt to quantify the error or 
accuracy of the results. This is a work in progress and we 
shall report in our future publications with a careful ex- 
amination of the different types of basis functions. And, 
calculate the dipole matrix elements in different gauges. 



1. One-valence 

Results from our calculations are listed in Table. IVII 
and component wise contributions are given in Table. 
One of the early works on the dipole matrix ele- 



VII 



ments of singly ionized alkaline-Earth metal ions is by 
Guet and Johnson [93 . The many-body method they 
used is MBPT and RPA, and numerical basis set. At 
the DF level the values of Guet and Johnson [93] are in 
good agreement, for both Sr + and Ba + , with our results, 
this is evident from the values listed in Table. IVIII Their 
work is the only one in the literature on the electric dipole 



matrix elements of Sr + and our results are higher. The 
difference could be largely attributed to the higher or- 
der core-polarization effects associated with the random- 
phase approximation (RPA). The RPA effects are incor- 
porated in the coupled-cluster but not to higher order as 
in an iterative RPA calculations. 

For Ba + there are several theoretical calculations of 
the electric dipole matrix elements. A careful study on 
the electric dipole transition is desirable as it a promising 
candidate for a novel parity non-conservation experiment 
|41) . In terms of the many-body method and single parti- 
cle basis set, the calculations of Sahoo and collaborators 
[P2] are closest to our approach. They estimate the upper 
bound on the error in the reduced dipole matrix element 
as 1.7%, which implies that our results have errors larger 
than this. 

For Yb + , the work of Safronova and Safronova [95] 
is the only previous study on the electric dipole matrix 
elements. Their calculations are based on the third order 
relativistic MBPT and the excellent matching between 
the length and velocity gauge results indicates the results 
are quite accurate. Our results are close to their results, 
however, at this stage we do not attempt to estimate the 
accuracy of our results. 



2. Two-valence 

Our results of the dipole matrix elements of Sr, Ba and 
Yb are listed in the tables Table. [VTTlp Cl There are very 
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TABLE V: Magnetic dipole HFS constant for the atomic systems Sr, Ba, and 1 Yb, using relativistic coupled-cluster 
theory. All values are in atomic units. 



State 






Coupled-cluster terms 




Other work 


Exp result 




DF 


P hfs -DF 


One-body H hfB 


Two-body H hfB 


Total value 






Atomic 87 Sr; [Kr]5s 2 
















5s5p 3 Pi 


-178.983 


-0.120 


-49.121 


0.002 


-228.222 


- 


-260.765(l) j 


5s5p 3 P 2 


-200.670 


0.106 


-47.045 


0.002 


-247.607 


— 


-212.085(5) j 


5s4d 3 £>i 


145.335 


0.098 


6.348 


0.001 


151.586 


— 


— 


5s4d 3 P 2 


-56.824 


0.076 


7.095 


0.001 


-49.654 


— 


— 


5s4d 3 D 3 


-133.930 


-0.040 


.194 


0.002 


-133.778 


— 


— 


5s4d 1 D 2 


17.441 


0.062 


9.643 


0.001 


27.145 


— 


— 


5s4d 1 P 1 


11.802 


-0.225 


4.366 


-0.002 


15.941 


— 


— 


Atomic 137 Ba; [Xe]6s 2 
















6s5d 3 Pi 


-588.432 


0.344 


-19.420 


-0.004 


-607.512 


-547 a 


-521 c 


6s5d 3 D 2 


397.451 


-0.588 


9.067 


-0.009 


405.921 


405 a 


416 c 


6s5d 3 P 3 


543.921 


0.189 


4.843 


-0.008 


548.945 


443 a 


457 c 


6s5d 1 D 2 


-148.545 


-0.459 


-54.967 


-0.006 


-203.977 


-102 a 


-82 d 


6s6p 3 Pi 


736.066 


-0.310 


221.943 


-0.004 


957.695 


1160 a 


1151 e 


6s6p 3 P 2 


806.032 


-0.152 


204.186 


-0.010 


1010.056 


845 a 


— 


6s6p 1 Pi 


—181.658 


—0.074 


-38.353 


0.009 


—220.094 


— 107 a 


— 109 


Atomic 173 Yb; [Xe]4/ 14 6s 2 
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-197.665 


0.004 


-906.806 
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6s6p 3 P 2 
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-0.292 
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-863.853 


-745 b 


-738 h 


6s5d 3 Di 


550.679 


0.035 


64.941 


0.002 


615.657 


596 b 


563(1)' 


6s5d 3 P 2 


-456.152 


-0.064 


-40.287 


0.003 
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-351 b 
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-454.431 


-0.027 


-22.482 
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-420 b 
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60 h 
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71.840 
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100(18)' 



a Ref erence [84] . b Reference [67]. 

c Reference[85 . d Reference 86 . 

c Reference [S7J. f Reference [55] . 

s Reference [55] . h Reference [50] . 

' Reference gj. j Reference^. 



few theoretical studies on the dipole matrix elements of 
Sr and those are not in very good agreement with ours. 
In the case of Ba, the previous theoretical calculations 
were done by Dzuba and Ginges [50] . Our results, listed 
IX are good agreement with ref. [3U] for the 
]0 1 1 6s 2 ^o) and (6s6p 3 Dj\ 



in Table. 

(6s6p 3 D 



\D\\6s5d d Dj 



However, there are large deviations for the matrix ele- 
ments involving the 6s6p Pi and 6s5d X P>2 states. For 
Yb, the (6s6p 1 Pi||I?||6s5d 3 I?2) is significantly different 
from the previous results. However, there is a large dif- 
ference between the previous results of Porsev and col- 
laborators [HE], an d Migdalek and Baylis [57] as well. 
Large relative deviations, compared to results of ref. [56"] . 
are also observed for the (6s6p 3 Pi||Z)||6s5<i 1 P>2) and 
(6s6p 3 P 2 ||P>||6s5d 1 P> 2 ). 



VIII. CONCLUSIONS 

In this paper we describe in detail the Fock-space rel- 
ativistic coupled-cluster method for the two- valence sys- 
tems. It is based on an all particle treatment and we 
demonstrate the excitation energies of Ba and Yb are 
on par with those of the previous relativistic CC cal- 
culations. The key point is, we have implemented the 
Fock-space CCT with an incomplete but quasi-complete 
model space comprising of the ns 2 , ns{n— l)d and nsnp 
configurations. This choice of model space is optimal 
to study the low-lying states of the alkaline-Earth metal 
atoms and other two- valence atoms like Yb and Hg. Most 
importantly, with this model space one may avoid diver- 
gences arising from the intruder states. This ought to 
be highlighted as the literature on iterative studies of 
two-valence systems is replete with accounts of intruder 
state induced divergences. We emphasize that there are 



19 



TABLE VI: Magnitude of the electric dipole transition am- 
plitude for 87 Sr+, 137 Ba+, and 173 Yb+ ions. 



Other works 

3.060° 
4.325" 
3.052" 
1.355" 
4.109" 



Ion 


Transition 


This work 


87 Sr+ 


5pi/2 - 


-> 5Si/2 


3.2180 




5f3/2 - 


-> 5Si /2 


4.9223 




5pi/2 — 


-J- 4d 3/2 






5^3/2 - 


-> 4d 3/2 


1.4217 




5^3/2 — 


4d 5 /2 


4.5942 


137 Ba+ 


6pi/2 - 


-> 6si/ 2 


3.1974 




6P3/2 - 


-> 6s 1/2 


5.0330 




6pi/ 2 — 


-> 5d 3/2 


3.0898 




6P3/2 — 


-> 5d 3/2 


1.2448 




6^3/2 — 


5d 5/2 


4.1347 


173 Yb + 


6pi/ 2 - 


-> 681/2 


2.9069 




6p 3 /2 - 


-> 6si/ 2 


4.5256 




6pi/ 2 — 


-> 5d 3/2 


3.6317 




6P3/2 — 


-J- 5d 3/2 


1.4918 




6P3/2 — 


5d 5/2 


4.8500 



3.009" 
1.312" 
4.057" 



3.300", 3.36(1)'', 
4.658", 4.73(3)'', 
3.11(3)'', 
1.34(2)*, 
4.02(7)'', 

2.731 d 
3.845 d 
3.782 d 
1.546 d 
4.769 d 



3.272 c 
4.614 c 
3.008 c 
1.313 c 
4.054 c 



"Reference |93] ■ 
'Reference T2l . 
c Reference l94l . 
Reference 95\ ■ 



few detailed relativistic many-body calculations of two- 
valence excitation energies and even less on properties 
calculations. Not surprisingly, among the published re- 
sults there is a wide variation of the many-body methods 
and basis sets used in the studies. Considering the grow- 
ing importance of alkaline-Earth metal atoms in precision 
experiments and possible applications, a detailed inves- 
tigations on the two-valence systems is timely. 

We have also developed a method based on CCT to 
compute properties from the CC wave functions of two- 
valence systems. This is perhaps an initial step towards 
systematic investigation of the structure and properties 
of two-valence systems with CCT, which has not been 
attempted before. Based on our scheme, the comptuta- 
tional cost of two-valence CC calculations is marginally 
higher than the one-valence calculations. And the addi- 
tional cost is in solving the S^ 2 ' cluster amplitude equa- 
tions. Total number of which is far less than the closed- 
shell and one-valence cluster amplitudes T and S^, re- 
spectively. So in terms of computational implementa- 
tions, there is no reason why two- valence CCT should not 
be the preferred method as in one- valence systems. The 
important and essential details of our schemes and impl- 
mentations are provided to highlight important physics 
issues in the two-valence Fock-space CCT. Breif descrip- 
tions of the method and extensive references are provided 
to aid interested researchers to implement CCT of two- 
valence systems. 

From the many-body theory perspective, CCT based 
structure and properties calculations is certainly an at- 



tractive choice. Prime reason being the topologically con- 
nected nature of the CC operators and exponential form 
of the wave operator ensures the condition of size exten- 
sivity. A basic requirement of a legitimate many-body 
theory. Further more, the non-perturbative character of 
the CC wave operator makes it an ideal choice. It must 
be emphasized that, all of these considerations are at the 
level of the many-body theory. However, the accuracy 
of the results also depends on other factors like single 
particle basis set considered. 

In the method we have developed, the one-valence 
coupled-cluster wave functions occur as an intermediate 
step. Using this we calculate the electric quadrupole HFS 
constant B of Sr+, Ba+ and Yb+. The HFS constant 
B of the 5p 3 / 2 2 P3/2 and 5d 3 / 2 2 D 3 / 2 states of 87 Sr+ 
and 137 Ba + , respectively are closer to the experimental 
data than the other theoretical results. For all the ions 
studied, the np 3 / 2 2 P 3 / 2 state has very large contribu- 
tions from the core polarization effects, which is part of 
(S^Htfs + c.c.) in the CC properties calculations. Simi- 
laraly, for bdj 2 Dj states of Yb + , there is a large contri- 
bution from the core-polarization effects. A careful ac- 
counting of the core-polarization effects is crucial for all 
the ions to obtain accurate values of B and this is particu- 
larly true for Yb + . We also calculate the magnetic dipole 
HFS constant A of Yb+, except for the 5d 5 / 2 state the 
results are in agreement with the other theoretical and 
experimental data. Similarly, we get reliable results of 
the dipole matrix elements of these ions. 

The results of the excitation energies of the two- valence 
sector calculations are in agreement with the previous 
CC results. For properties calculations with CC wave 
functions, our results show deviations from the previous 
works. However, it must be mentioned that there have 
been very few attempts at theoretical properties calcula- 
tions of two- valence systems. And the previous works are 
based on MCDF or using a collage of single particle wave 
functions. The later may require finer analysis for pre- 
cision studies as the linked-cluster theorem, which forms 
the basis of many-body theory, is based on a uniform sep- 
aration of the total Hamiltonian into zeroth order and 
perturbation. This is not the case when orbitals with 
different central potentials are used in the calculations. 
With MCDF method, a large scale structure and proper- 
ties calculations of neutral atoms with high Z is plagued 
with convergence issues. Among all the states, the sin- 
glet states nsnp 1 P\ and ns(n — l)d l D 2 require further 
attention as the properties involving these state exhibit 
largest deviations from other theoretical results. A simi- 
lar pattern is observed in the other theoretical results as 
well. Based on our studies and careful analysis, the ob- 
served deviations of the two-valence properties may be 
attributed to the basis, V N ~ 2 potential, we have used. 
We expect calculations with V N ~ 1 potential basis could 
improve the results. 

In conclusion, relativistic Fock-space coupled-cluster 
theory has theoretical and computational advantages for 
structure and properties calculations of two-valence sys- 
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terns. In this article we report the development of an 
all particle two-valence relativistic Fock-space coupled- 
cluster theory and have demonstrated a scheme for prop- 
erties calculations with the CC wave functions. 
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TABLE VII: The electric dipole transition amplitude, contributions from different terms in the coupled-cluster theory. 
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TABLE VIII: El transition amplitudes for the atomic system 


87 Sr, using relativistic coupled-cluster theory. All values are in 
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-0.0001 


0.8257 


0.00001 




0.4497 


0.16°, 0.162 6 


x Pi^ 


'So 


-4.2442 


0.0000 


0.5283 


-0.00002 




-3.7159 


5.28", 5.238\ 


1.9539 c 


3 P ^ 


3 D 1 


2.6323 


-0.0002 


-0.3131 


0.00000 




2.3190 






3 Px^ 


3 D 1 


2.2652 


0.0013 


0.1116 


0.00000 




2.3781 






3 P 2 ^ 


3 D 1 


0.5849 


0.0009 


0.3772 


-0.00001 




0.9630 






x Pi^ 


3 D 1 


0.2255 


-0.0005 


-0.2496 


-0.00001 




-0.0247 






3 A^ 


3 D 2 


3.9538 


0.0012 


-0.5437 


-0.00001 




3.4114 






3 P 2 ^ 


3 D 2 


2.2646 


0.0006 


0.0700 


0.00001 




2.3352 






x Pi^ 


3 D 2 


0.7531 


-0.0004 


-0.1054 


0.00001 




0.6473 






3 P 2 ^ 


3 D 3 


-5.3938 


-0.0001 


0.3111 


0.00001 




-5.0828 






3 P X ^ 


'D 2 


-0.8822 


0.0012 


0.0387 


0.00000 




-0.8423 


0.19" 




3 P 2 -^ 




-0.2854 


0.0002 


-0.3597 


-0.00001 




-0.6448 


0.10" 




x Pi^ 


X D2 


-4.5484 


-0.0008 


0.3860 


-0.00001 




-4.1632 


1.92" 





"Reference 
'Reference 
c Reference 
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TABLE IX: El transition amplitudes for the atomic system Ba, using relativistic coupled-cluster theory. All values are in 
atomic units. 



Transition Coupled-cluster terms Other work 







DF 




One-body Hhis 


Two-body Hhf B 


Total value 




3 P1- 




0.3888 


-0.0003 


-0.8090 


0.00000 


0.4205 


0.4537" 


X P1- 


-> 


-4.6768 


-0.0002 


0.6730 


0.00000 


4.0040 


5.236 a 


3 P - 


-> 3 P! 


2.6203 


0.0004 


-0.2585 


0.00000 


2.3622 


2.3121 a 


3 P1- 


-> 3 /A 


2.2405 


-0.0023 


-0.0195 


0.00000 


2.2187 


2.0108" 


3 D 

P 2 — 


v 3 n 


—0.5715 


0.0019 


n /(too 

—0.4782 


0.00000 


1.0478 


0.5275 




-> 3 D, 


-0.3364 


-0.0002 


0.3546 


0.00000 


0.0180 


0.1047 a 


3 Pl- 


-> 3 P 2 


3.8886 


-0.0022 


-0.1585 


0.00001 


3.7279 


3.4425 a 


3 P 2 ~ 


-> 3 £ 2 


-2.2265 


0.0013 


0.1007 


-0.00001 


2.1245 


2.024° 


X Pl- 


-> 3 P 2 


-0.3874 


-0.0005 


0.0833 


-0.00000 


0.3046 


0.4827 a 


3 P 2 ~ 


"> 3 £>3 


5.3410 


-0.0004 


-0.3409 


-0.00001 


4.9997 


4.777 a 


3 Pl- 


-> X P 2 


-1.1039 


-0.0018 


0.1178 


0.00000 


0.9879 


0.1610 a 


3 P 2 ~ 


-> X P 2 


0.4458 


0.0005 


0.5219 


0.00001 


0.9682 


0.1573" 


J Pl- 


-»■ ^ 


4.4933 


-0.0011 


-0.0773 


0.00000 


4.4149 


1.047" 



"ReferencelgDl. 



TABLE X: El transition amplitudes for the atomic system 173 Yb, using relativistic coupled-cluster theory. All values are in 
atomic units. 

Transition Coupled-cluster terms Other work 







DF 


Hhfs-DF 


One-body Hhis 


Two-body Hms 


Total value 




3 P1- 




0.1445 


-0.0003 


-0.5320 


0.00000 


-0.3878 


0.54(8)", 0.44\ 0.587 c 


X P1- 


-> % 


-3.8641 


-0.0001 


0.5999 


-0.00001 


-3.2643 


4.40(80)", 4.44 6 , 4.89 d , 4.825 c 


3 Po- 


-> 3 D 1 


2.7296 


0.0001 


-0.3209 


0.00000 


2.4088 


2.61(10)", 2.911 c 


3 P1- 


-> 3 D 1 


2.3473 


-0.0005 


0.1811 


0.00003 


2.5279 


2.26(10)" 


3 P 2 - 


-> 3 D 1 


-0.5997 


0.0000 


-0.2343 


0.000002 


-0.8340 


0.60(12)" 


X P1- 


-> 3 D 1 


-0.4503 


-0.0002 


0.1702 


0.00000 


-0.2803 


0.27(10)", 0.24 6 


3 P1- 


-> 3 D 2 


3.9875 


-0.0005 


-0.6480 


0.00002 


3.3390 


4.03(16)" 


3 P 2 - 


-> 3 £ 2 


-2.2940 


-0.0002 


-0.1010 


-0.00003 


-2.3952 


2.39(1)" 


X P1- 


-> 3 P 2 


0.0716 


-0.0003 


0.0660 


0.00000 


0.1373 


0.32(6)", 0.60 6 


3 P 2 - 


-> 3 £ 3 


5.6130 


0.0000 


-0.3215 


-0.00001 


5.2915 


6.12(30)" 


3 P1- 


-> X P 2 


-1.1920 


-0.0006 


0.0995 


0.00000 


-1.0931 


0.54(10)" 


3 P 2 - 


-> ^ 2 


0.5946 


0.0002 


0.3601 


0.00000 


0.9549 


0.38(8)" 


X P1- 


-> X P 2 


4.7006 


-0.0002 


-0.5209 


0.00000 


4.1795 


3.60(70)" 



"Reference |96] . 

6 Reference E ■ 
"Reference 68] . 

Reference |98 . 



